Removal of a single photon by adaptive absorption 
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We present a method to remove, using only linear optics, exactly one photon from a field-mode. 
This is achieved by putting the system in contact with an absorbing environment which is under 
continuous monitoring. A feedback mechanism then decouples the system from the environment as 
soon as the first photon is absorbed. We propose a possible scheme to implement this process and 
provide the theoretical tools to describe it. 



I. INTRODUCTION 



In the last decades experimental technology has 
reached the point in which quantum effects cannot only 
be observed but also controlled. Sensors have become 
very precise and individual quantum systems can be ma- 
nipulated through coupling to semiclassical driving fields. 
Despite this, we are still far from having full control over 
quantum systems. In the field of quantum optics, the two 
main reasons for this are that the coupling to the driv- 
ing fields is in practice quite limited, and the non-desired 
coupling to the environment produces damping and de- 
coherence effects which hinder controllability. Classical 
control theory has been quite successful in manipulating 
classical systems with complex or partially unknown dy- 
namics. Classical control consists of closely monitoring 
the system, processing the measured data and applying 
the appropriate feedback signal to modify the evolution 
of the system. In the field of quantum optics this idea 
has been used to stabilize the phase or intensity of lasers. 
However, this reduction of noise could be explained clas- 
sically. 

Non-linear effects are used to reduce noise in the quan- 
tum domain. Yamamoto et al [|l| presented the first 
scheme which relied only on feedback to generate pho- 
tocurrents with noise below the standard quantum limit. 
Based on Langevin equations, Shapiro et al. Q gave the 
first quantum and semiclassical descriptions of feedback 
control. Wiseman and coworkers || derived a closed-loop 
master equation describing the instantaneous feedback of 
a homodyne current measured outside an optical cavity. 
This work, based on quantum trajectories, made clear the 
connection between control and quantum measurement 
theory and made possible a complete understanding of 
'quantum-limited' feedback. Subsequently several appli- 
cations and elaborations of this first quantum theory of 
continuous feedback appeared. In general the feedback 
does not have to be instantaneous and can be a functional 



of the entire history of measurement results, thus typi- 
cally resulting in a non-Markovian evolution. Doherty 
and Jacobs Q showed that in determining the form of 
this functional it is useful to formulate feedback control 
in two steps: estimation of the state of the system, and 
application of appropriate control inputs to affect the dy- 
namics. In H Doherty et al. extended this idea further 
to investigate how methods in the well developed classi- 
cal control theory can apply quantum feedback control. 
From the application point of view, feedback has been 
proposed as a means to improve atom localization in 
standing waves J(||| , cooling of cavity mirrors j?J , con- 
trolling coherence of two-level atoms fla] , inverting quan- 
tum jumps [0], line narrowing in atomic fluorescence [fiol , 
slowing down of decohcrcncc effects jnj, and preserving 
macroscopic superpositions of light fields . 

In most of these applications, feedback is used to keep 
an open system under control. The back-action of mea- 
surement in the 'quantum-limited' feedback is nearly con- 
sidered as a drawback, and the feedback is enforced by 
the action of a driving field that shapes the potential. In 
our proposal the feedback itself is used to create an other- 
wise difficult operation. In fact, the evolution is tailored 
only by the continuous measurement procedure. 

The idea in this work is to show that feedback control 
implemented by simple linear optical elements can be of 
great use in manipulating quantum states of light. This 
is, at present, particularly interesting in the context of 
quantum communication where photons are the only se- 
rious qubit carriers. Their very weak coupling strengths 
suppress decoherence effects, but at the same time ren- 
der rather difficult the implementation of quantum gates. 
Any tools which perform non-trivial tasks without rely- 
ing on photon-photon interactions are therefore very de- 
sirable. The use of linear optical elements has already 
proven to facilitate some tasks Here, we address 

the particular problem of removing precisely one photon 
from a field mode. Our proposal is based on a very simple 
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form of feedback and for its implementation only linear 
optical elements are needed. The central idea is that of 
adaptive absorption and consists of coupling the mode to 
an environment and monitoring that environment for the 
absorption of a photon. Once a photon has been detected 
the coupling to the environment is switched off. If the 
coupling is maintained until a photon has been absorbed, 
then this procedure removes exactly one photon from the 
field. 

The paper is organized as follows. Section || briefly 
reviews the main properties of linear absorption. The 
master equation describing the evolution of the field in 
the absorber is introduced. These results are used in 



Section II] where we give the evolution of the field in 
adaptive absorption, and in Section |y| where we give 
a possible implementation. In Section |V], we study the 
evolution under adaptive absorption for some relevant 
input states. Section VI treats the question obtaining 
information about the input field from the time at which 
the absorption occurs. Finally, in Section VII we present 



some further elaborations and conclude the paper. 



II. LINEAR ABSORPTION 

Absorption in a linear medium is characterized by a 
reduction of optical field amplitude by a factor 772 which 
is independent of the field amplitude. A coherent state 
\a) after being partially absorbed during its propagation 
through, for example, an optical fiber is transformed to 
rjia^. This is formally the same relation as between 

the input and one of the output arms of a beam splitter. 
The P-representation (see Eq. (17) below) allows us to ex- 
tend the previous formal analogy to all input states. The 
beam-splitter also serves as a model for phase-sensitive 
damping [jl4| , p~5| by feeding the second input port with a 
squeezed state. 

The action of beam splitters is widely studied (see for 
example |l(| ) . Some properties which are relevant in our 
context are: a) The output modes are in general entan- 
gled and therefore the transmitted state is a mixed state, 
b) Coherent states are an exception to the previous state- 
ment, c) When the vacuum and n photons are fed into 
the input ports, the photon number distribution in the 
output ports follows the same statistical law for distribut- 
ing n classical particles: i.e. from n photons the absorber 
will let through m photons at random with probability r\ 
for each, while absorbing (m — n) with probability 1 — r;. 
The resulting probability distribution for the transmitted 
field is, 



(1) 



d) Vacuum fluctuations from the input port are added to 
the the transmitted field. This illustrates the fact that 



dissipation is always accompanied by extra fluctuations. 

The beam splitter analogy gives us a relation between 
the input and output photonic states. While the field is 
in contact with the absorbing medium, its evolution is 
non-unitary because of the losses, but may be described 
by the master equation 



^ = T(2ap{t)J 



a) ap{t) — p(t)a'a 



(2) 



where a is the annihilation operator of the field-mode. 
This is the master equation of a harmonic oscillator cou- 
pled to a broad-bandwidth bath ]l7| . It is useful to sep- 
arate the evolution superoperator into two terms defined 
by 



J p = apa' 

hp = — — \o)ap + pa) a] 



(3) 
(4) 



The master equation (j^) can now be rewritten in the 
compact form, 



p = 2T(J + L)p. 



(5) 



A solution of this differential equation may be expressed 
in a number of ways including p8[ , 



p(t) =J2 Mm I dt rn-l ■ • • / dtie 2 ^'"'-) 



2rJe 2r£( * m -*'— l) • • • 2rje 2r£tl /3(0) 



(6) 



This formal solution allows for an interpretation in 
terms of conditioned evolution in agreement with quan- 
tum measurement theory. The solution is formed from 
terms of the form, 



p c (t) = exp[2rL(i - t N )}2TJeyip[2TL(t N 
•••2rjexp[2r J Lii] j o(0). 



t N -i)}- 



(7) 



Each occurrence of J corresponds to the loss of a pho- 
ton. The exponential terms account for the non-unitary 
evolution of the system between photon absorptions. It 
follows that 



Pc(t) 



Pc(t) 

Tr(pc(t)) 



(8) 



can be interpreted as the evolved density matrix condi- 
tioned by the loss of N photons to the environment at 
times ti, ti, . . . , 1 7v . The norm of p c {t) gives the proba- 
bility of this particular sequence of detections to occur. 
The unconditional evolution described in Eq. is then 
obtained by averaging over all possible histories, i.e. sum- 
ming over the number of lost photons m, and integrating 
over their corresponding absorption times t\ . . .t m . The 
action of the superoperator J is sometimes referred to as 
a 'quantum jump', since it produces a finite change in 
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the state in an infinitesimal time. In contrast, when no 
quantum jump occurs, the change is infinitesimal but not 
unitary. 

A coherent state after the loss of a photon and after a 
time At without losses is transformed, respectively, as, 



IV. 



J \a)(a\ 
3 \a){a\ 



|a) (a\ 



— g-M 2 (l-exp(-2rAt)) 

x |ae- rA *Yae- rA *| 



(9) 
(10) 



One immediately notices that for an initial coherent 
state p(0) = |a)(a|, each integrand in Eq. (||) is propor- 
tional to |ae _r *)(ae _r * | . Hence, the unconditioned state 
at a time t is also |ae _rt ). This shows the equivalence of 
the master equation (Q) to the beam splitter from above 
(with rj = e _2r * .) The master equation description and 
its interpretation in terms of conditioned evolution gives 
us a tool to formalize adaptive absorption in the following 
section. 



III. ADAPTIVE ABSORPTION 

In this section we show how to remove a single pho- 
ton from an field-mode prepared in any state. For this 
we need to couple our field to the broad-bandwidth envi- 
ronment as in the previous section. The state will then 
evolve according to Eq. (^|), or Eq. (Q) if we are able 
to keep track of the jump times. Now suppose that we 
switch off the coupling to the environment immediately 
after the a single photon has been absorbed (see below for 
a possible implementation). In this case the conditional 
state after a photodetection at time t\ < t will be 



p c (t) = 2TJe^ Ltl p(0) 



(11) 



We can predict the evolution of this system in which 
a photon may be absorbed at any time before i, or no 
photon is absorbed. The corresponding unconditioned 
normalized density matrix is, 



p(t) = exp[2rXi]p(0) 



dti2rjexp[2rXii]p(0). (12) 



It is worth noting that this evolution is explicitly non- 
Markovian. The future evolution depends not only on 
the state but also on its previous history. This can also 
be seen by differentiating the previous equation, 



p(t) = 2T(J + L)e 2rLt p{0) 

= 2T(J + L)p n0 detection 



(13) 



where the 'no detection subscript denotes the density op- 
erator for which no photon has been absorbed. 

Before examining the properties of Eq. (|l2]), we de- 
scribe a scheme to implement the adaptive absorption. 



IMPLEMENTATION OF ADAPTIVE 
ABSORPTION 



Vacuum 



INPUT 





N 

FIG. 1. Beam splitter model for a field-mode coupled to an 
environment 

The basic idea is depicted in Fig. [l]. The input signal is 
send through a series of weak beam splitters. Their reflec- 
tion coefficient is very small, so that there is a negligible 
probability for more than a single photon to be reflected. 
A photodetector is placed in the weakly coupled arm of 
each beam splitter. What we have described until now, is 
our field mode coupled to the broad-bandwidth environ- 
ment. The corresponding evolution is then described by 
Eq. (^|) or Eq. (0) depending on whether the firing times 
of the detectors are known or not. In the latter case the 
input-output relation is equivalent to a single beam split- 
ter of finite reflectivity. In particular this means that all 
vacuum contributions interfere to give a single vacuum 
contribution, so that the act of continuously monitoring 
does not entail extra fluctuations. 

We seek to enforce the conditional evolution. This 
means that if a signal is registered in the first detector 
Di, no further beam splitters should be interposed. On 
the other hand, if the first detector does not fire, the light 
is sent to the second beam splitter and the photodetector 
Du is checked for a firing signal. The process is repeated 
until a photodetector fires. 

In order to implement this procedure, the experimental 
arrangement in Fig. || might be used. 
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, k INPUT 



FIG. 2. Possible experimental setup for single photon ex- 
traction by adaptive absorption. 
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In this scheme only four optical elements and a single 
photodetector are needed. When a light pulse is trapped 
in the loop enclosed between the modulated beam split- 
ter MBS, the mirrors Mi and M2, and the flip mirror, 
it is effectively coupled to the environment through the 
modulated beam splitter MBS. The value of its reflec- 
tivity is set close to one and can be slightly increased via 
an electric current to make it fully reflective. In order 
to effectively feed the initial state into the loop the flip 
mirror has to open. 

The photodetector keeps track of the photon losses and 
as soon as it fires, it sends a feedback signal to the flip- 
mirror so that the pulse can exit the loop without any 
further losses. Alternatively, the coupling of the field 
to the environment can be switched off by turning the 
MBS fully reflective. The flip mirror is then switched 
whenever the output is desired. This last alternative al- 
lows the system to provide the unconditioned evolved 
states. Otherwise different sub-ensembles will emerge at 
the output at different times. Depending on the response 
times of the photodetector, control electronics, modu- 
lated beam-splitter, and the flip-mirror one can add a 
delay line within the loop, for example between the two 
mirrors Mi and M%. 

Any real experiment will inevitably be only an approx- 
imation to the model described above. In practice every 
beam splitter has a finite reflectivity and has some in- 
ternal losses, and detectors have efficiencies below one. 
Inclusion of such effects will degrade the performance of 
the device. 

V. PROPERTIES OF THE EVOLVED STATE 



time until the first jump occurs. The conditioned state 
at the output when the detector fires at t\ is given by 

p«(*) = |ae- rtl X«e" rtl |- ( 14 ) 
This occurs with the probability 

p{h\a)At = 2r|a| 2 e" 2rtl exp[-|a| 2 (l - e~ 2rtl )}dt. (15) 

Following Eq. (|l2|) the unconditional evolution of the 
coherent state is 

p a (t) = exp[-|a| 2 (l - e- 2r *)] |ae- rt )(ae- rt | 

+ f dii2r|a| 2 e- 2rtl 
Jo 

x exp[-|a| 2 (l - e" 2 " 1 )] \ae~ rtl )(ae~ rtl | . (16) 

The P-representation of this state, defined by the equa- 
tion 

P a (t) = J d 2 [3P a (P,t)\P){P\, (17) 

is given by 

P Q (/3,i) = {exp[-|a| 2 (l-e- 2r *)]|/3|- 1 
x S(\(3\ - \a\e- Ft ) + 2e^ 2 -^H(\/3\ - He~ r *) 
x H(\a\-\p\)}8(axga-aig0), (18) 

where H (x) is the Heaviside step function. The existence 
and form of this P-function allow us to draw some con- 
clusions: 1) The state produced in this way will have 
essentially classical properties. 2) The removal of the 
photon does not change the argument of the coherent 
state in that the P-function is zero for phases other than 
arg a. 



In this section we list some properties of the evolved 
state. Eqs. ([ll]) and ( [l2|) give the conditional and uncon- 
ditional evolution for an arbitrary initial state and we 
now present the evolution for two important examples. 



A. Coherent state 

We consider first the important example of an initial 
coherent state p Q (0) = As shown in Eq. (||), 

a coherent state remains unaffected under the action of 
the jump superoperator J. This can be understood in the 
context of the implementation in Fig. ||. As already men- 
tioned, coherent states have the peculiarity that when 
sent through a beam splitter the two outcoming beams 
are not entangled. This implies that the transmitted 
beam is the same independently of the measurement out- 
come at the photodetector. In any case, while coupled 
to the environment, the coherent amplitude is damped 
following Eq. (|l0|). Thus the output field is determined 
only by the time which the states spend in contact with 
the environment (inside the loop in Fig. |^), that is the 



P a (P,t) 




2e 2n lal 2 lal 2 P 

FIG. 3. P-distribution for a coherent state |a) evolved until 
time t (full line). P°"(f3,t) is non- vanishing only at the fixed 
phase arg (3 — arg a. 

The plot of the P-function P a ((3,t) in Fig. || gives us 
an intuitive picture of the evolution of the unconditioned 
state. The initial state \a) transforms after a time t into 
a mixture of coherent states with damped amplitudes 
\(3\ G [|a|c~ rt , a|). The contribution in the mixture of 
each coherent state also decays exponentially. The singu- 
lar peak (* in Fig. |J) corresponds to the case where the 
detector has not fired, and it takes the required weight 
to keep P a {(3, t) normalized. 
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B. Number states and photon statistics 

Another important basis-state is the number state \n). 
With probability 



p(t 1 \n)dt= 2Tne- 2Tntl dt 



(19) 



the number state \n) will lose one photon at time t\, to 
produce the conditioned state 



p n c {t) = \n-l){n-l\ 



(20) 



for times t > t\. This state is independent of the time 
t\ in which the jump occurs. As a result, the uncondi- 
tioned state when no information on the detection events 
is available, reduces to the simple form, 



P n (t) 



-2nTt 



\n)(n\ + (1 



-2nTt 



)|n-lXn-l|. (21) 



In contrast to the coherent state case, the number state 
soon evolves into a pure state. The probability of not 
extracting one photon from \n) falls exponentially with 
time, thus the extraction is bound to occur given enough 
time. After long times t ~ -» oo a single photon is removed 
from the signal deterministically (unless, of course, there 
was no photon to begin with). 

This last remarks are in fact valid for any initial state 
p(0) = Y^ n '=o Pn,n> \ n )( n '\- Indeed, the photon number 
probability distribution at a time t is given by 



Pn{t) = (n\p(t) \n) 



-2nTt 



P»(0) + (1 



-2nTt 



Wl(O). 



And for t — > oo we have 



p„(oo) = p n +i(0) + Pa(0)S n ,o 



(22) 



(23) 



which means that the probability distribution is shifted 
one step to smaller values of n. The initial vacuum com- 
ponent is of course not shifted and still receives the con- 
tribution from the initial one-photon component. In this 
limit of large times the mean and variance of the photon 
number can be easily calculated using Eq. (p3|), 

fi(oo)=n(0)+po(0)-l, (24) 
An 2 (oo) = Ati 2 (0) -po(0)[2n(0) +Po(0) - 1] 

= An 2 (0) -po(0)[n(0) +n(oo)]. (25) 

The mean photon number decreases by an amount 
equal to the initial probability of having a non-vacuum 
signal. The variance decreases by an amount dictated 
by the vacuum probability and also by the mean pho- 
ton number. From the results in the previous section 
we know that adaptive absorption cannot produce non- 
classical light from initially classical light, in the sense 
that it preserves the positivity of the P-function. De- 
spite this, it still enables us to change substantially the 



photon statistics of the incoming field. The normally or- 
dered variance of the number operator is 

: An 2 : (oo) = Aro 2 (oo) — n(oo) 

= :An 2 : (0) + 1 - po(0)[2n(0) + po(0)]. (26) 

It is now easy to find an example in which an initial 
super-Poissonian distribution (: An 2 :> 0) gives a sub- 
Poissonian output field. An incoming field p A (0) with 
Pn(0) = Po^ n ,o + (1 — Po)$n,N has a normal ordered vari- 
ance : An 2 := N(l — po)(Npo — 1) which is positive for 
N > Pq 1 ■ After long times the corresponding output 
field will have the normal ordered variance, 



:An 



2 .A 



(oo) = (JV-l)(l-po)[(JV-l)p -l] 1 (27) 



which becomes negative for 1 < N < p$ + 1 . Summariz- 
ing, for every non-integer value of p^ 1 there exists an inte- 
ger N which fulfills the inequalities p$ 1 < N < p$ 1 + 1 so 
that the corresponding super-Poissonian input field be- 
comes sub-Poissonian after the adaptive absorption pro- 
cedure. 



VI. ADAPTIVE ABSORPTION AS PHOTON 
NUMBER MEASUREMENT 

Along the course of this paper we have investigated 
adaptive absorption as a way to manipulate photonic 
states. By now we have expressed in many ways the 
idea that the induced evolution is non-unitary and that 
the coupling to the environment plays a crucial role in 
this manipulation. In quantum information theory this 
type of evolution is always accompanied by an outflow 
of information to the environment. In this section we 
show that indeed adaptive absorption can also be un- 
derstood as a procedure to measure the input photonic 
state. The balance between information gain and distur- 
bance is highly non-trivial and delicate. Moreover, when 
the system under study is an optical field-mode we find 
practical limitations in that many of the interesting quan- 
tum operations require non-linear optical devices. Here 
we present a potentially valuable tool which only makes 
use of linear optical elements. 



In Section VB we have seen that for long times, the 
effect of adaptive absorption on the photon number dis- 
tribution was to shift it down by one. Now we will see 
that by keeping track of the detection times t\ we can 
obtain information on the initial photon-number. Taking 
into account that the probability of removing one photon 
from a state of n photons is the same as for the analogous 
classical problem, it is easy to understand that the larger 
the initial photon number, the easier (i.e., the faster) it 
is to extract a single photon. The time t\ at which the 
photon is absorbed is our measurement outcome. In fact, 
at each time step we are effectively performing the mea- 
surement defined by the positive operator-valued measure 
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or POVM (H) {ni = 2dtr a t a ,n = 1 - Hi}. Since the 
absorption POVM Hi is innnitcsimally small, most of 
the times we get the "0" outcome. The value of t\ con- 
veys the result of the measurements up to t\. In order to 
know the information obtained from the knowledge of t\ , 
the central quantity is the probability p(n\t a ) of having 
a certain number n of photons in the input given that 
the measurement outcome is t\ — t a . By invoking Bayes 
theorem we can write this conditional probability as 



P(n\t a ) 



p{t a \n)p n 

P(ta) 



p(t a \n) 



Em=lP(*«l n )' 



(28) 



where p(t a ) — E^Lo PnP(t a \n) and we assume for defi- 
niteness a flat a priori probability distribution p n = p 
Vn. Use of the conditional probability given in Eq. (ITT 
gives 



P(n\t a ) 



-2T{n+l)t a r 2Tt a 



1 



(29) 



We note that summation of this expression over n from 
1 to oo gives unity as it should. The exponential depen- 
dence on n may provide the means to discriminate be- 
tween states with widely differing photon numbers with- 
out detecting all of the photons (see Figure ||) . The tech- 
nique cannot, of course, determine the number of photons 
without error. The benefit of this technique as a measure- 
ment, however, is likely to be in those situations where 
we can only afford to weakly perturb the mode. A possi- 
ble application of this idea to quantum cryptography is 
discussed briefly in the concluding section. 




0.5 1 
FIG. 4. p{n\t a ) for V 
n = {l,2,5}. 
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the possibility of removing more than one photon. The 
fact that one can stop the absorption process as soon 
as one photon is absorbed is the novel aspect of this 
work. For example this possibility has given rise to a 
novel eavesdropping attack on quantum key distribution 
(QKD) p0[ |. There are two factors that might render 
insecure many of the current implementations of QKD 
pl[ : the signals are weak coherent pulses instead of sin- 
gle photons and the communication channels have very 
high losses. The eavesdropper can acquire the secret key 
whenever he succeeds in extracting a single photon from 
the multiphoton part of the signal. Adaptive absorption 
allows the eavesdropper to take advantage of this security 
gap by using current technology. 

By repeatedly applying the single-photon extraction 
one can extract any number of photons in a controlled 
way. Moreover, as we have seen in the previous sec- 
tion, in process of extracting single-photons we also ob- 
tain information on the input system. Thus adaptive 
absorption implements weak measurements, where the 
measured system is only slightly disturbed. 

Adaptive absorption is a simple example of a general 
idea that can give rise to many other applications. An- 
other example, that falls outside the scope of this paper, 
is adaptive amplification as a means to add a single ex- 
citation into a mode. This general idea is to use the re- 
sults of a continuous measurement to modify the system 
Hamiltonian or its coupling to the measurement appa- 
ratus. This line of action may lead to novel quantum 
operations or ease the realization of operations with too 
high technological demands. 
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VII. CONCLUSIONS 

In this paper we have presented the idea of extracting 
a single photon from an arbitrary single-mode field via 
adaptive absorption. The extraction of a single photon 
can also be achieved by post-selection, that is by setting 
up a beam splitter and considering at the end only those 
situations in which a single photon was removed. Of 
course this is not deterministic and would in fact lead to 
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